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PREFACE 


This  report  describes  the  results  of  a  detailed  examination  of 
the  radiation  and  heat  budgets  of  the  Rand  version  of  the  Mintz-Arakawa 
two-level  atmospheric  general  circulation  model  for  a  January  control 
run.  The  Mintz-Arakawa  model  is  used  extensively  in  the  Rand  Climate 
Dynamics  Project,  for  the  Defense  Advanced  Research  Projects  Agency, 
and  this  study  represents  part  of  the  continuing  effort  to  verify  and 
improve  the  model.  Considerations  of  the  sort  undertaken  here  are 
necessary  to  insure  that  the  radiation  and  heat  terms  in  the  model  are 
good  approximations  to  the  real  atmosphere,  and  that  we  may  reasonably 
perform  numerical  experiments  with  the  model.  In  addition,  this  work 
provides  us  with  reference  levels  for  the  analysis  of  such  experiments. 

(A  similar  report  dealing  with  the  July  radiation  and  heat  budget  is 
planned. ) 

Reports  related  to  this  study  include  the  following:  C.  Schutz 
and  W.  L.  Gates,  Global  Climatic  Data  for  Surface ,  800  mb,  400  mb: 
January ,  R-915-ARPA,  November  1971;  W.  L.  Gates,  E.  S.  Batten,  A.  B. 
Kahle,  and  A.  B.  Nelson,  A  Documentation  of  the  Mintz-Arakawa  Two-Level 
Atmospheric  General  Circulation  Model ,  R-877-ARPA,  December  1971;  C. 
Schutz  and  W.  L.  Gates,  Supplemental  Global  Climatic  Data:  January , 
R-915/1-ARPA,  May  1972;  W.  L.  Gates,  The  January  Global  Climate  Simu¬ 
lated  by  the  Two-Level  Mintz-Arakawa  Model:  A  Comparison  with  Obser- 
vatton,  R-1005-ARPA,  November  1972;  W.  L.  Gates,  The  January  and  July 
Climates  Simulated  by  a  Global  2-Level  General  Circulation  Model:  A 
New  Comparison  with  Observation ,  1973  (in  preparation);  Staff,  Climate 
Dynamics  Project,  The  Frvtd  General  Circulation  Model ,  1973  (in 
preparation). 


SUMMARY 


The  radiative  terms  and  heat  budget  in  the  Rand  version  of  the 
Mintz-Arakawa  two-layer  atmospheric  general  circulation  model  are 
examined  and  compared  with  the  available  values  from  observational 
and  theoretical  studies.  The  absorption,  reflection,  and  scattering 
experienced  by  the  solar  radiation  in  the  atmosphere  and  at  the  earth's 
surface,  the  long-wave  radiative  fluxes  at  the  top  and  bottom  of  the 
model  atmosphere,  the  nonradiative  transfer  of  heat  via  sensible  heat 
and  evaporation  and  condensation,  and  the  planetary  albedo  are  dis¬ 
cussed  and  utilized  in  deriving  the  heat  budgets  across  the  boundaries 
of  the  model  and  for  the  atmospheric  column  and  the  earth-atmosphere 
system.  In  addition,  the  horizontal  transports  of  heat  necessary  to 
balance  these  budgets  are  derived.  These  considerations  are  all 
limited  to  zonal  and  global  results  for  Januan  from  a  control  inte¬ 
gration  of  the  model. 

The  model  indicates  that  too  much  solar  radiation  reaches  and  iB 
absorbed  at  the  surface  due  primarily  to  insufficient  absorption  in 
the  atmosphere  and,  in  a  few  areas,  insufficient  reflection.  As  for 
the  long-wave  fluxes,  in  the  tropical  latitudes  the  flux  divergences 
over  the  atmospheric  column  (the  difference  between  the  net  flux  at  the 
surface  and  the  net  flux  at  the  top  of  the  atmosphere)  appear  to  be  too 
large.  In  the  mid-latitudes  of  both  hemispheres  the  values  are  in 
better  agreement  with  comparable  values,  while  at  higher  latitudes  they 
again  diverge.  It  should  be  noted,  however,  that  since  both  the  flux 
of  solar  radiation  reaching  the  surface  and  the  long-wave  fluxes  are 
strongly  dependent  on  the  distribution  of  cloudiness,  moisture,  and 
temperature,  which  still  need  improvement,  it  is  not  evident  that  the 
fault  lies  necessarily  in  the  radiative  portions  of  the  model. 

The  various  components  of  the  heat-budget  calculations  seem  rea¬ 
sonable  when  compared  with  those  of  other  investigators.  The  globally 
averaged  values  of  the  horizontal  heat  transport  do  not  go  to  zero,  how¬ 
ever,  and  certain  adjustments  to  the  model  are  suggested.  The  Mintz- 
Arakawa  heat-budget  components  are  generally  the  largest  discussed, 
suggesting  that  the  model  possesses  a  rather  vigorous  general  circulation. 

Preceding  page  blank 


In  general,  the  various  comparisons  between  the  results  of  the 
Mintz-Arakawa  model  and  those  of  other  studies  indicate  that  the  model 
is  realistic.  Certain  modifications  are  suggested  as  being  immediately 
practical,  while  other  improvements  are  seen  to  depend  on  improvement 
of  the  primary  meteorological  parameters  upon  which  the  radiation  cal¬ 
culations  are  based. 
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I.  INTRODUCTION 


This  report  describes  a  detailed  examination  of  the  radiative 
terms  and  heat  budget  in  the  Rand  version  of  the  Mintz-Arakawa  two- 
layer  atmospheric  general  circulation  model.  Such  a  study  was  under¬ 
taken  in  order  to  identify  what  areas  of  the  radiation  modeling  need 
improvement.  In  addition,  a  preliminary  attempt  has  been  made  to  ascer¬ 
tain  why  these  areas  are  weak.  As  part  of  Rand's  program  in  climate 
dynamics,  this  radiation  study  will  be  used  as  a  reference  or  control 
in  the  analysis  of  future  experiments  on  the  mechanisms  of  climate 
change . 

In  this  report  the  radiative-energv  terms  as  calculated  from  con¬ 
trol  experiment  14  for  January  will  be  described.  In  addition,  so  that 
the  entire  heat  budget  may  be  examined,  the  latent  and  sensible  heat- 
flux  terms  will  also  be  considered.  In  this  version  of  the  model  the 
criteria  for  determining  cloudiness  have  been  changed,  January  rather 
than  mean  annual  sea-surface  temperatures  used,  and  several  minor 
errors  removed  (Gates  et  al.  ,  1971;  Staff,  Climate  Dynamics  Project, 
1973).  We  will  compare  the  results  with  a  representative  sample  of 
observational  and  theoretical  studies.  Included  will  be  studies  which 
are  primarily  concerned  with  the  radiative-energy  budgets  and  therefore 
utilize  climatological  values  for  the  relevant  meteorological  param¬ 
eters,  and  also  general  circulation  models  which  generate  their  own 
values  for  these  parameters  as  well  as  the  radiative  quantities.  This 
work  complements  the  studies  of  the  observed  January  climate  (Schutz 
and  Gates ,  1971,  1972),  the  simulated  climate  for  January  using  an 
earlier  version  of  the  model  (Gates,  1972),  and  an  extensive  review  of 
the  January  climate  simulated  by  this  newer  version  of  the  control  ex¬ 
periment  (Gates,  1973). 

As  used  here,  January  refers  to  the  30-day  period  from  December  31 
to  January  29.  Except  in  a  few  cases  where  data  were  saved  only  every 
six  hours,  the  January  averages  are  based  on  numbers  calculated  every 
half  hour.  In  particular,  all  the  energy  terms  involved  in  the  heat- 
budget  calculations  used  in  this  report  are  based  on  half-hour  values. 


Since  some  of  the  studies  available  for  comparison  only  dealt  with 
seasonal  averages,  we  have  found  it  necessary  in  those  cases  to  com¬ 
pare  our  January  values  with  average  values  for  the  period  from  Dec¬ 
ember  through  February.  We  have  confined  our  study  to  zonal  and  global 
averages.  (A  subsequent  report  will  deal  with  the  July  and  annual 
radiation  budgets.) 

In  order  to  proceed  with  this  investigation  of  the  January  radia¬ 
tion  and  heat  budget  produced  by  the  Mintz-Arakawa  model,  it  is  neces¬ 
sary  to  look  first  at  the  distributions  of  cloudiness,  water  vapor, 
temperature,  evaporation,  precipitation,  sensible  heat,  and  surface 
albedo  as  generated  or  used  by  the  model.  The  examination  of  these 
parameters  is  only  of  secondary  importance  to  this  study;  as  already 
mentioned,  a  more  extensive  review  of  these  parameters  from  a  per¬ 
spective  not  limited  Just  to  the  radiative  and  heat  budgets  is  under 
way  and  will  be  reported  on  elsewhere  by  Gates.  Tne  distributions  of 
these  quantities,  which  are  present  here  for  the  purpose  of  comparison, 
will  generally  be  referenced  to  the  radiation  studies  in  which  they 
were  used  rather  than  to  their  original  sources.  This  initial  discus¬ 
sion  is  followed  by  an  examination  of  the  solar  and  long-wave  radiation 
terms  available  from  the  model.  Finally,  the  radiative  and  nonradia- 
tive  components  of  the  heat  budget  are  considered.  Zonally  averaged 
values  of  all  these  quantities,  as  derived  from  the  Mintz-Arakawa 
model,  are  given  in  the  Appendix. 


II.  METEOROLOGY 


One  of  the  most  important  meteorological  parameters  in  the  model 
as  far  as  radiation  is  concerned  is  the  cloudiness.  In  an  earlier 
version  of  the  model  the  cloudiness  was  much  too  low,  permitting  too 
much  solar  radiation  to  reach  and  be  absorbed  by  the  surface  of  the 
earth  (Gates,  1972)  and  too  much  long-wave  radiation  to  escape  to 
space.  The  criteria  for  determining  cloudiness  have  been  revised  by 
Koenig  (see  Staff,  Climate  Dynamics  Project,  1973)  and  the  cloud  amounts 
are  now  significantly  improved.  The  zonal  average  of  the  total  cloud¬ 
iness  from  the  Mintz-Arakawa  model  is  shown  in  Fig.  1  (see  also  Gates, 
1973).  Values  of  cloudiness  given  oy  London  (1957),  Katayama  (1967), 
van  Loon  (19/2),  and  the  data  both  of  ETAC  (1971)  and  of  Miller  (1970) 
as  presented  by  Schutz  and  Gates  (1971),  are  also  shown,  London's 
and  Katayama 's  climatological  values  were  prepared  for  their  radiation 
studies,  and  the  others  as  parts  of  climate  studies.  Haurwitz  (1972) 
and  Sasamori  et  al.  (1972)  both  used  van  Loon's  cloud  data  for  their 
radiative  studies.  While  the  Mintz-Arakawa  cloudiness  still  appears 
to  be  too  small  on  a  global  basis,  the  low  values  in  the  tropics  corre¬ 
spond  rather  well  with  the  recent  Tiros  satellite  data  of  Miller,  which 
show  considerably  less  cloud  cover  in  those  latitudes  than  do  most  of 
the  previous  estimates.  In  the  higher  northern  latitudes  the  various 
distributions  of  cloudiness  vary  rather  consistently  with  latitude, 
with  the  Mintz-Arakawa  model  having  slightly  less  cover  than  most  of 
the  others.  The  agreement  between  the  data  sets  in  the  southern  hemi¬ 
sphere  is  much  poorer.  This  is  particularly  true  south  of  50°S,  where 
the  Mintz-Arakawa  cover  is  smallest  and  each  distribution  shows  a  dif¬ 
ferent  trend  with  increasing  latitude.  The  effect  of  cloud  distribution 
on  the  radiation  xs  quite  noticeable  and  will  be  discussed  below. 

The  radiation  is  also  affected  by  the  water  present  in  the  atmos¬ 
phere.  Zonally  averaged  precipitable  water  amounts,  compiled  by  Starr 
et  al.  (1957)  as  part  of  a  climatological  study,  aro  given  in  Fig.  2. 

The  distributions  used  by  the  radiative  studies  of  London  (1957) , 


In  all  the  figures  the  values  generated  by  the  general  circulation 
nodels  are  denoted  by  a  double  asterisk,  while  values  generated  by  the 
radiation  models  are  denoted  by  a  single  asterisk.  Climatologically 
derived  values  are  unmarked. 


Latitude  (deg) 

Fig.  1-The  zonally  averaged  total 


averaged  precipftable  water  vapor. 


Katayama  0-967),  and  Sasamori  et  al.  0972)  are  also  presented  in 
Fig.  2.  Haurwitz  0972)  used  essentially  the  same  distribution  as 
that  by  Sasamori  et  al.  Rather  than  generating  their  own  values, 

Holloway  and  Manabe  0971)  used  Telegadas  and  London's  (1954)  values 
for  the  lower  troposphere  and  Murgatroyd's  (1960)  for  the  upper  tropo¬ 
sphere.  In  the  tropics  the  values  generated  by  the  Mintz-Arakawa  model 
differ  substantially  from  the  other  values,  which  are  all  in  relatively 
close  agreement.  The  hydrological  cycle  is,  indeed,  felt  to  be  one  of 
the  weakest  areas  of  the  model,  and  work  is  under  way  to  improve  it. 

The  temperature  of  the  atmosphere  is  important  in  both  the  long¬ 
wave  radiative  calculations  and  the  energy  conversions  involving  mois¬ 
ture.  Figure  3  compares  the  January  zonally  averaged  global  ground 
temperatures  generated  by  the  Mintz-Arakawa  model  with  the  climato¬ 
logical  surface-air  values  collected  by  Schutz  and  Gates  (1971)  from 
Crutcher  and  Meserve  (1970)  for  the  northern  hemisphere  and  from 
Taljaard  et  al.  (1969)  for  the  southern  hemisphere.  The  latter  values 
were  used  by  Sasamori  et  al.  (1972)  and  Haurwitz  (1972)  in  their  studies. 
Also  included  are  the  climatological  surface-air  values  prepared  by 
London  (1957)  for  his  heat -balance  study  of  the  northern  hemisphere. 

The  distributions  are  similar  in  the  tropical  and  middle  latitudes. 
However  at  higher  latitudes  the  model  generates  excessively  high  tem¬ 
peratures.  In  the  Antarctic  this  is  probably  due  to  an  albedo  error. 
While  this  was  corrected  prior  to  the  January  simulation,  residual 
effects  are  still  noted. 

The  surface  albedo  of  the  model  is  shown  in  Fig.  4,  along  with  the 
albedos  given  by  Katayama  (1967)  and  Sasamori  et  al.  (1972).  The 
albedo  used  in  the  Mintz-Arakawa  model  is  based  on  the  mean  surface 
albedo  for  January  of  Posey  and  Clapp  (1964)  as  tabulated  by  Schutz 
and  Gates  (1971),  while  the  albedo  of  Sasamori  et  al.  is  based  on  the 
work  of  Landsberg  et  al.  (1965)  and  Budyko  (1963).  There  seems  to  be 
fairly  good  agreement  among  the  various  data. 

The  zonally  averaged  evaporation  rates  for  January,  as  derived  by 
the  general  circulation  models  of  Mintz-Arakawa  and  Holloway  and  Manabe 
(1971)  and  as  collected  by  Budyko  (1963)  and  presented  by  Schutz  and 
Gates  (1971),  are  given  in  Fig.  5.  Figure  6  presents  the  precipitation 


Latitude  (deg) 

Fig.  3-The  zonal ly  averaged  ground  and  surface  air  temperature. 


ly  averaged  surface  albedo. 


Fig.  5-The  zonally  averaged  evaporation. 
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Fig.  6-The  zonally  averaged  precipitation. 


-11- 


rates  for  the  two  genera],  circulation  models,  as  well  as  the  climato¬ 
logical  data  of  Moller  (1951)  as  given  by  Schutz  and  Gates  (1972).  ^or 
both  rates  the  Mintz-Arakawa  model  shows  almost  twice  as  much  moisture 
conversion  in  the  tropics  as  do  the  other  distributions.  It  is  evi¬ 
dent  that  the  hydrologic  cycle  in  the  Mintz-Arakawa  model  proceeds  at 
much  too  rapid  a  pace. 

In  order  to  complete  the  consideration  of  the  heat  budget,  we  must 
include  the  sensible  heat  exchange.  The  flow  of  sensible  heat  may  occur 
either  to  ot  from  the  surface  of  the  earth.  When  we  examine  the  zonally 
averaged  distribution  of  sensible  heat  flux  in  Fig.  7,  as  calculated  by 
the  Mintz-Arakawa  model,  it  appears  that  the  flux  is  downward  over  as 
much  of  the  earth  as  it  is  upward.  Budyko  (1963),  as  collected  by 
Schutz  and  Gates  (1971),  and  Holloway  and  Manabe  (1971),  however,  show 
only  an  upward  transfer  of  sensible  heat  in  their  zonal  distributions. 


Fig.  7-The  zonally  uveraged  sensible  heat  flux. 
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III.  SHORT-WAVE  RADIATION 


The  incoming  solar  radiation  used  by  the  model  for  January  is  based 
°Y  solar  constant  of  2.0  ly  min"-  (cal  cnf2  min"1),  and  calculated  as 
a  i. unction  of  the  time-dependent  solar  orbital  parameters  and  local  solar 
zenith  angle.  Since  the  total  solar  radiation  reaching  the  top  of  the 
atmosphere  is  unaffected  by  the  model,  it  should  closely  resemble  other 
distributions  of  incident  solar  radiation,  being  dependent  only  on  the 
value  chosen  for  the  solar  constant  and  whether  January  or  the  three 
winter  months  (December,  January,  and  February)  are  being  charactered. 
The  zonal  averages  of  the  January  incident  radiation  specified  by  the 
Mintz-Arakava  model  are  shown  in  Fig.  8.  Also  shown  for  the  northern 
hemisphere  are  values  from  Katayama  (1966,  1967)  for  January,  based  on 
a  solar  constant  of  1.94  ly  min"1,  and  from  London  (1957),  based  on  a 
solar  constant  of  2.0  ly  min"1  and  averaged  over  the  three  winter 
months.  For  the  southern  hemisphere  the  values  of  Sasamori  et  al. 

r 


(1972)  for  January  15,  based  on  a  solar  constant  of  2.0  ly  min 
giver . 
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are 


Of  the  incoming  solar  radiation,  roughly  15  percent  is  absorbed 
by  the  atmosphere  (which  is  thus  directly  heated) ,  35  percent  is  re¬ 
flected  back  to  space  by  the  clouds,  the  atmosphere  (back  scattering), 
and  the  ground,  and  50  percent  is  absorbed  by  the  surface  of  the  earth 
The  latter  energy  may  be  stored  temporarily  in  the  form  of  an  increase 
in  the  surface  temperature  but  is  eventually  lost  by  the  earth's  sur¬ 
face  in  the  form  of  long-wave  radiation,  sensible  heat,  and  latent 
heat,  thus  heating  the  atmosphere  indirectly. 

In  the  Mintz-Arakawa  model,  atmospheric  absorption  of  solar  radia¬ 
tion  is  calculated  for  each  layer  successively  as  a  function  of  the 
water-vapor  content  of  that  layer.  Clouds  are  incorporated  simply  by 
assigning  them  as  equivalent  water-vapor  content.  All  other  gaseous 
and  aerosol  absorption  is  neglected.  On  the  other  hand  London  (1957) 
reports  on  the  atmospheric  absorption  by  ozone  and  clouds  individually 
and  water  vapor  and  dust  together.  Katayama  (1967)  gives  the  absorp¬ 
tion  due  to  water  vapor,  dust,  and  clouds  separately,  and  Sasamori 
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et  al.  (1972)  treat  ozone,  water  vapor,  and  clouds  separately  and 
molecular  oxygen  and  carbon  dioxide  together.  The  results  of  London 
and  Sasamori  et  al.  are  given  for  the  entire  atmosphere,  and  those  of 
Katayama  for  the  troposphere.  Thus  their  results  cannot  be  compared 
with  the  Mintz-Arakawa  values  calculated  for  the  individual  layers. 

The  zonally  awraged  values  of  atmospheric  absorption  for  the 
various  data  sets  are  given  in  Fig.  9.  It  should  be  noted  that  because 
the  absorption  of  solar  radiation  is  one  of  the  smallest  components  of 
the  radiation  budget,  the  scale  of  this  figure  is  five  times  larger 
than  that  used  in  the  rest  of  the  figures  in  this  report.  The  data 
from  Kataya-  are  for  the  troposphere  only.  The  upper  curves  for 
Mintz-Arakawa,  London,  and  Sasamori  et  al.  represent  the  total  absorp¬ 
tion  in  the  atmosphere.  The  lower  curve  for  Mintz-Arakawa  is  for  the 
troposphere  only,  the  difference  between  the  two  curves  being  due  to 
the  small  amount  of  absorption  by  water  vapor  in  the  stratosphere. 

The  lower  curves  for  the  data  of  London  and  Sasamori  et  al.  do  not 
include  ozone  absorption.  Thus,  if  we  assume  that  all  ozone  absorp¬ 
tion  occurs  in  the  stratosphere  and  all  other  absorption  in  the  tropo¬ 
sphere,  the  lower  curves  represent  the  troposheric  absorption  curves 
for  London  and  Sasamori  et  al. ,  to  be  compared  with  the  data  of  Katayama 
and  the  lower  Mintz-Arakawa  curve. 

If  we  first  look  at  this  "tropospheric"  absorption,  it  is  apparent 
that  the  Mintz-Arakawa  model  absorbs  considerably  less  solar  radiation 
in  the  lower  and  middle  latitudes  than  is  indicated  by  the  other  in¬ 
vestigators.  This  is  particularly  true  when  the  results  are  compared 
with  those  of  Sasamori  et  al.  for  the  southern  hemisphere  and  i3  at 
least  in  part  due  to  the  low  amount  of  water  vapor  in  the  Mintz-Arakawa 
model.  However,  referring  to  Fig.  2,  it  appears  that  the  Mintz-Arakawa 
water-vapor  amounts  are  low  between  20°N  and  30°S,  while  the  tropo¬ 
spheric  absorption  is  low  between  30°N  and  60°S.  This  suggests  that 
more  fundamental  deficiencies  in  the  model  may  exist.  One  may  be  the 
lack  of  absorbers  other  than  water  vapor.  Another  may  be  the  assump¬ 
tion  that  only  that  part  of  the  solar  radiation  of  wavelength  greater 
than  0.9  p  is  subject  to  absorption.  Others  (see,  for  instance, 
Dopplick,  1970)  have  divided  the  spectrum  at  0.7  y. 


le  zonally  averaged  solar  radiation  absorbed  by  the  total  atmosphere  (up 
id  the  troposphere  (lower  curves). 
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The  tropospheric  absorption  ls  the  quantity  of  interest  uhen 
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a  sorption  of  solar  radiation  In  the  Mlntz-Arakawa  model  means  that  the 

-unt  of  solar  radiation  incident  on  the  cloud  tops  is  increled  in 

comparison  with  the  amounts  present  In  the  other  studies,  fhe  solar 

tadiat  on  in  the  Mintz-Arakawa  model  has  traversed  one-quarter  of  the 

2d  ;rSPhere  bef"e  rea‘hi"8  the  9  percent  of  the 

c  ends,  one-half  the  mass  before  reaching  the  next  44  percent  of  the 

C  »u  a.  and  three-quarters  of  the  mass  by  the  time  It  reaches  the  final 
percent  of  the  clouds.  Thus,  not  only  the  lack  of  D?n  k 
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6  reduced  tropospheric  absorption  will 

This  t  T  7  enhanCed  8°lar  radiation  reaching  the  cloud  tops, 
his  t  nds  to  compensate  for  the  reduced  cloud  cover,  so  that  the  re¬ 
acted  radiation  is  not  too  different  from  that  found  by  others.  In 


the  southern  hemisphere  in  particular,  starting  with  this  premise  that 
more  radiation  is  reaching  the  cloud  tops  in  the  Mintz-Arakawa  model 
than  in  the  other  studies,  the  distributions  of  reflected  solar  radia¬ 
tion  correlate  well  with  the  corresponding  cloud  distributions.  Around 
35*S  and  near  the  south  pole  the  Mintz-Arakawa  model  has  about  the  same 
or  slightly  more  cloud  cover  than  van  Loon,  and  the  reflected  flux  is 
greater  than  that  of  Sasamori  et  al.  In  the  tropics  the  Mintz-Arakawa 
model  has  considerably  less  cloud  cover  but  only  slightly  smaller 
amounts  of  reflected  radiation.  Only  between  ^5°S  and  75°S,  where  the 
difference  between  the  smaller  Mintz-Arakawa  cloud  cover  and  that  of 
Sasamori  et  al.  is  greatest  and  the  absorption  discrepancy  is  somewhat 
reduced,  is  there  a  significantly  smaller  reflected  flux  in  the  Mintz- 
Arakawa  model. 

The  zonally  averaged  distributions  of  solar  radiation  incident 
upon  the  earth's  surface  in  January  are  shown  in  Fig.  11.  In  addition 
to  the  values  produced  by  the  Mintz-Arakawa  model,  the  results  of  London, 
Sasamori  et  al.,  Katayama,  and  Budyko  (1963)  as  given  by  Schutz  and 
Gates  (1971)  are  also  presented.  As  would  be  expected,  given  the  lower 
amounts  of  absorption  found  for  the  Mintz-Arakawa  model,  the  model's 
solar  radiation  incident  upon  the  surface  is  larger  than  any  of  the 
values  reported  by  the  other  studies.  In  middle  and  high  latitudes  of 
the  winter  hemisphere  the  various  distributions  are  quite  similar,  not 
so  much  because  greater  accuracy  has  been  achieved  there  but  rather 
simply  because  the  hemisphere  is  relatively  dark.  The  details  of  the 
curves  in  Fig.  11  can  be  explained  quite  well  on  the  basis  of  our  pre¬ 
vious  observations  about  the  solar  radiation  in  the  Mintz-Arakawa  model. 
In  the  tropics  there  is  insufficient  absorption  in  the  model,  and  there¬ 
fore  too  much  radiation  reaching  the  ground.  From  30°S  to  35°S  the 
solar  radiation  reaching  the  ground  is  about  the  same  in  the  model  as 
in  the  other  data  sets.  In  this  region  the  model  reflects  more  than 
the  others  due  to  a  fairly  high  cloud  cover,  and  this  high  reflection 
compensates  for  the  insufficient  absorption.  Proceeding  towards  higher 
southern  latitudes,  the  absorbed  amounts  remain  smaller  in  the  model 
while  the  relative  amount  of  reflected  radiation  decreases,  so  the 
curves  in  Fig.  11  diverge  with  more  radiation  reaching  the  ground  in 
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t;he  model.  The  largest  discrepancies  are  reached  around  60°S  to  70°S, 
where  the  Mintz-Arakawa  cloudiness  reaches  a  relative  minimum  and  both 
absorption  and  reflection  are  low. 

Figure  12  presents  the  zonally  averaged  January  values  of  that  part 
of  the  incident  solar  radiation  absorbed  by  the  earth's  surface.  The 
same  studies  are  reported  on  as  in  the  previous  figure  except  that  Budyko 
values  are  replaced  by  those  of  the  Geophysical  Fluid  Dynamics  Labora¬ 
tory  (GFDL)  model  (Holloway  and  Manabe,  1971).  A  comparison  of  Figs. 

11  and  12  clearly  illustrates  the  strong  effect  the  large  snow-surface 
albedo  values  have  in  the  high  latitudes  of  the  southern  hemisphere. 

The  planetary  albedo  is  defined  as  the  ratio  of  the  solar  radia- 
t  ’-eflected  and  scattered  to  space  to  the  solar  radiation  incident 
on  the  earth-atmosphere  system.  Zonally  averaged  values  for  January  as 
calculated  by  Katayama  (1967),  Sasamori  et  al.  (1972),  and  the  Mintz- 
Arakawa  model  are  given  in  Fig.  13.  In  addition,  recent  satellite  mea¬ 
surements  by  Vonder  Haar  (1972)  and  Raschke  et  al.  (1973)  are  included. 

As  pointed  out  by  Vonder  Haar  and  Suomi  (1971),  the  satellite  observa¬ 
tions  tend  to  yield  lower  albedos  than  the  other  studies  do.  The 

Planetary  albedo  of  the  Mintz-Arakawa  model  appears  to  lie  within  the 
range  of  the  other  values. 

In  summary,  then,  comparing  the  Mintz-Arakawa  model  with  other 
data,  we  find  insufficient  absorption  and,  in  a  limited  region,  insuf¬ 
ficient  reflection  of  the  solar  radiation  as  it  passes  through  the 
atmosphere.  This  allows  too  much  solar  radiation  to  reach  and  be  ab¬ 
sorbed  at  the  earth's  surface.  When  modifying  the  model,  care  should 
be  taken  to  correct  the  cloud,  water-vapor,  and  ozone  amounts,  and 
possibly  also  the  albedo  values  simultaneously. 


zonally  averaged  planetary  albedo. 
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IV*  LONG-WAVE  RADIATION 


values  of  the  net  outgoing  lo„g-„ave  radiation  flux  at  the  200  mb 
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Although  the  Mintz-Arakawa  distribution  of  the  net  surface  flux 
has  a  form  somewhat  like  that  of  the  other  profiles  presented  in  Fig. 

15,  it  differs  substantially  from  them  in  various  ways.  For  thejriost 
part  its  values  are  largest,  being  anywhere  from  10  to  30  ly  day 
larger  than  the  other  values.  In  the  equatorial  latitudes,  where  the 
differences  are  greatest,  the  amount  of  precipitable  water  generated 
by  the  Mintz-Arakawa  model  is  considerably  lower  than  the  climatologi- 
cally  derived  amounts  used  in  all  the  other  studies.  The  fact  that 
the  Mintz-Arakawa  atmosphere  has  less  moisture  than  the  others  means 
that  the  downward  flux  at  the  surface  will  be  smaller  and  accordingly 
the  net  upward  flux  greater.  At  higher  latitudes  the  agreement  among 
the  water-vapor  concentrations  is  much  better,  as  is  the  agreement 

between  surface  flux  distributions. 

The  cloud  cover  generated  by  the  Mintz-Arakawa  model  in  the  tropi¬ 
cal  latitudes  is  approximately  20  percent  lower  than  the  climatological 
values  used  by  the  other  studies.  Also  it  is  approximately  constant 
between  20°N  and  20°S,  while  the  climatological  values  show  a  slight 
peak  around  the  equator.  Thes,  differences  will  certainly  contribute 
to  the  discrepancies  between  the  Mintz-Arakawa  net  surface  flux  di¬ 
stribution  and  those  calculated  by  others. 

It  is  also  possible  that  differences  in  the  surface  temperatures 
may  account  for  some  of  the  discrepancies  between  the  net  surface  flux 
profile  of  Mintz-Arakawa  and  those  of  the  other  studies.  The  surface 
temperature  profiles  available  are  given  in  Fig.  3.  Only  in  the  polar 
latitudes  do  the  differences  appear  to  be  significant.  There  the  Mintz 
Arakawa  model  produces  temperatures  as  much  as  30  degrees  warmer  than 
climatological  values.  The  existence  of  these  warmer  polar  surfaces 
will  contribute  to  the  larger  net  outgoing  fluxes  found  by  the  Mintz- 
Arakawa  model.  In  addition  the  model  produces  unrealistically  large 
diurnal  temperature  variations,  which  are  not  apparent  in  the  30-day 
averages.  Since  the  outgoing  radiative  flux  has  a  T  dependence, 
larger  values  of  average  outgoing  flux  may  be  generated  even  though 

the  average  ground  temperatures  are  the  same. 

In  the  southern  hemisphere  south  of  the  polar  front  the  Mintz- 
Arakawa  model  first  shows  an  increase  in  the  net  flux  and  then  a 
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decrease  beyond  75°S.  This  pattern  is  clearly  tied  to  the  distribution 
of  cloudiness  in  the  high  latitudes,  where  the  cloud  cover  first  de¬ 
creases  and  then  increases.  The  same  feature  exists  in  the  northern 
hemisphere  but  is  much  less  well  developed  because  there  is  less  cloud 

cove...  None  of  the  other  studies  show  such  extreme  variations  of  the 
net  long-wave  radiation. 

In  the  polar  latitudes  of  the  southern  hemisphere  the  net  surface 
fluxes,  as  calculated  by  Sasamori  et  al.  (1972)  and  Haurwitz  (1972), 
are  only  20  ly  day"1  lower  than  their  peak  hemispheric  values.  This 
quite  substantial  increase  with  increasing  latitude  is  a  function  of 
the  decrease  in  cloud  cover  and  water-vapor  content  toward  the  pole. 
Holloway  and  Manabe  (1971)  show  a  smaller  recovery,  and  the  Mintz- 

Arakawa  values,  as  already  noted,  have  a  local  maximum  not  at  the  pole 
but  at  approximately  75° S. 

Figure  16  presents  the  zonally  averaged  values  of  the  net  infrared 
flux  divergence  of  the  atmospheric  column  as  calculated  by  the  Mintz- 
Arakawa  model  and  by  the  other  investigations  discussed  in  association 
with  Fig.  15.  These  values  are  of  interest  because  they  provide  a  mea¬ 
sure  of  the  net  infrared  cooling  in  the  column.  It  should  be  pointed 

out  that  the  atmospheric  column  does  not  extend  to  the  same  level  in 
all  cases. 

Generally  the  lower  layers  of  the  atmosphere  contain  enough  water 
vapor  and  carbon  dioxide  to  be  relatively  opaque  in  the  infrared  part 
of  the  spectrum.  Thus  since  the  spread  between  the  effective  emitting 
temperatures  of  the  upward  and  downward  fluxes  passing  through  the 
bottom  of  the  column  is  not  large,  the  net  flux  at  the  surface  will  be 
relatively  small.  For  this  reason  the  meridional  variations  of  the  net 
flux  divergence  are  primarily  a  function  of  the  flux  at  the  top  of  the 
column.  We  have  already  noted  that  the  net  flux  at  the  top  of  the 
atmosphere  decreases  with  increasing  latitude  because  of  the  corre¬ 
sponding  decrease  in  temperature.  The  same  explanation  accounts  for 
the  decrease  of  the  net  flux  divergence  with  increasing  latitude. 

In  the  tropical  latitudes  the  Mintz-Arakawa  results  are  largest, 
in  part  at  least  because  of  the  small  cloud  cover  generated  in  those’ 
latitudes.  One  might  have  anticipated  that  Holloway  and  Manabe' s  (1971) 


averaged  long-wave  radiative  flux  divergence  for  various 


-31- 


values  should  have  exceeded  those  of  our  model,  since  their  column 
extends  to  the  top  of  the  atmosphere.  However  since  they  used  London's 
(1957)  climatological  values  of  cloud  cover,  which  are  nearly  20  per¬ 
cent  larger  than  those  of  the  Mintz-Arakawa  model  in  the  tropics,  this 
is  not  the  case.  At  higher  latitudes,  where  the  cloud  cover  is  in 
somewhat  better  agreement,  Holloway  and  Manabe's  values  are  indeed 
greater  than  those  of  the  Mintz-Arakawa  model.  Actually,  in  the  middle 
and  high  latitudes  of  the  northern  hemisphere  the  agreement  between 
London  (1957),  Katayama  (1967),  and  Mintz-Arakawa  is  quite  good. 

In  the  southern  hemisphere  there  is  a  greater  disparity  between 
the  various  distributions.  The  values  of  Sasamori  et  al.  (1972)  are 
probably  higher  than  they  should  be  as  the  result  of  their  using  a 
normalization  procedure  which  implicitly  assumes  that  there  is  no  over¬ 
lapping  between  clouds  at  different  levels.  Unfortunately  the  cloud 
values  used  by  Haurwitz  (1972),  while  permitting  the  clouds  to  overlap, 
are  based  on  an  as  yet  poorly  developed  climatological  record,  so  it 
is  difficult  to  estimate  the  validity  of  those  results  either.  Finally, 
since  Holloway  and  Manabe's  (1971)  cloud  cover  was  taken  from  season¬ 
ally  transposed  values  of  London  (1957) ,  those  results  are  equally 


ambiguous . 

Both  the  GFDL  and  Mintz-Arakawa  general  circulation  models  produce 
distributions  of  n^fc  flux  divergence  centered  not  at  the  equator  but 
at  approximately  1(TS  as  would  be  expected  during  the  southern  hemi¬ 
sphere  summer.  In  addition  the  net  flux  divergences  of  Sasamori  et  al. 


(1972)  for  the  southern  hemisphere  are  fairly  constant  from  about  25°S 
to  the  equator,  If  the  cloud  distributions  used  by  these  three  studies 
are  examined  (Fig.  1),  it  can  be  seen  that  all  have  lower  values  between 


10°S  and  15°S  than  at  the  equator,  which  would  yield  the  sort  of  dis¬ 
tribution  of  net  flux  divergence  we  find.  However,  a  similar  dip 
appears  at  10°N  to  20*N  in  both  of  those  cloud-cover  distributions 
which  extend  into  the  northern  hemisphere,  and  this  suggests  that  there 
should  be  a  relative  minimum  in  the  flux  divergence  at  the  equator, 
something  we  do  not  find  in  any  of  the  distributions.  Warmer  temper¬ 
atures  are  skewed  slightly  toward  the  southern  hemisphere  and  undoubtedly 
contribute  to  this  distribution  of  the  net  flux  divergence  in  the  tropics. 
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In  conclusion  it  would  appear  that  in  the  tropical  latitudes  of 
the  Mintz-Arakawa  model  too  much  long-wave  radiation  is  escaping  through 
the  tropopause.  If  the  entire  tropospheric  column  is  considered,  this 
excessive  loss  is  in  part  compensated  for  by  the  overestimation  of  the 
net  upward  flux  at  the  surface.  However,  the  flux  divergence  over  the 
entire  column  is  still  somewhat  too  large.  At  higher  latitudes  in  the 
wrthern  hemisphere  the  flux  divergences  look  reasonable,  although  near 
the  pole  the  surface  temperatures  are  probably  too  warm.  In  the  mid¬ 
latitudes  of  the  southern  hemisphere  the  values  ar_  acceptable,  but  at 
higher  latitudes  the  model  does  not  seem  to  function  as  well.  It  should 
be  pointed  out,  however,  that  because  of  the  strong  dependence  of  the 
long-wave  fluxes  on  the  distributions  of  clouds,  moisture,  and  temper¬ 
ature,  to  criticize  the  radiative  sections  of  the  model  would  be 
premature. 
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V.  RADIATION  AND  HEAT  BUDGET 


In  this  section  the  net  radiative  fluxes  through  the  upper  (tropo- 
pause)  and  lower  surfaces  of  the  Mintz-Arakava  model,  as  well  as  the 
nonradiative  heat-flux  terms  at  the  lower  surface,  will  be  examined. 

The  heat  budgets  for  the  model  atmospheric  column  and  earth-atmosphere 
system  will  also  be  considered.  From  the  calculated  heat  budget  for 
the  atmospheric  column  the  amount  of  heat  that  has  to  be  transported 
horizontally  to  balance  that  budget  will  be  derived.  In  addition, 
despite  the  fact  that  the  ocean  temperature  is  fixed,  values  for  the 
ocean  transport  of  heat  required  by  the  Mintz-Arakawa  model  will  be 
calculated  under  the  assumption  that  the  values  given  by  Newell  et  al. 
(1969)  for  the  heat  storage  in  the  ocean  can  be  used.  The  balances 
and  their  various  components  will  also  be  compared  with  similar  quan¬ 
tities  given  by  others. 

The  zonal  distribution  of  the  January  flux  through  the  top  of  the 
model  atmosphere  (200  mb),  shown  in  Fig.  17,  is  comprised  of  the  in¬ 
coming  solar  radiation  minus  the  total  reflected  solar  radiation  and 
the  net  outgoing  long-wave  radiation.  Although  this  flux  is  taken  at 
the  tropopause,  it  will  be  quite  similar  to  the  flux  at  the  top  of  the 
atmosphere.  Even  though  it  is  not  necessarily  true  for  any  given  lat¬ 
itude,  London  (1957)  shows  that  globally  the  radiative  fluxes  at  the 
tropopause  and  the  top  of  the  atmosphere  are  equal,  with  the  ozone 
absorption  of  solar  radiation  just  balancing  the  stratospheric  net 
loss  of  long-wave  radiation  to  space.  The  net  flux  at  the  tropopause 
is  cl imato logically  important,  since  the  geographical  distribution  of 
this  energy  imbalance  is  the  ultimate  driving  force  for  the  atmospheric 
and  oceanic  circulations. 

In  addition  to  the  flux  at  the  tropopause  from  the  Mintz-Arakawa 

m°del*  Flg'  17  8lves  the  flu“s  in  thn  northern  hemisphere  at  the  top 
of  the  atmosphere  from  London  (1957)  and  at  200  mb  from  Katayama  (1967) 
and  In  the  southern  hemisphere  at  the  top  of  the  atmosphere  from 
Sasamori  et  al.  (1972),  as  well  as  the  global  satellite  data  of  Vender 
Hear  and  Suomi  (1971).  while  following  the  general  trend  of  the  other 
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distributions,  the  Mintz-Arakawa  values  are  slightly  smaller  throughout 
most  of  the  low  latitudes  and  quite  a  bit  larger  between  50°S  and  75°S. 

In  the  low  latitudes,  while  both  the  model's  net  outgoing  long-wave 
and  incoming  short-wave  radiation  are  larger  than  the  values  reported 
by  the  studies  with  which  they  are  compared,  the  differences  are  greater 
in  the  case  of  the  long-wave  radiation.  Accordingly,  the  values  of  the 
model's  net  downward  flux  are  somewhat  smaller  than  those  of  the  other 
studies  presented  in  Fig.  17.  The  large  peak  in  the  high  southern 
latitudes  occurs  in  a  region  where  the  cloud  cover  and  albedo  of  the 
model  are  relatively  low  compared  with  the  other  values  available. 

Thus  the  amount  of  solar  radiation  reflected  to  space  (see  Fig.  10)  is 
also  relatively  low.  In  addition,  the  net  outgoing  long-wave  flux  at 
200  mb  is  comparatively  small  in  that  region.  The  cumulative  effect 
is  to  produce  an  excessively  large  value  for  the  net  flux  entering  the 
troposphere  in  these  latitudes.  This  will  undoubtedly  affect  the 
atmospheric  circulation  in  the  high  southern  latitudes,  a  region  where 
it  is  felt  that  the  model's  performance  is  less  satisfactory  (see,  for 
Instance,  Gates,  1972). 

We  next  look  at  the  heat  budget  across  the  lower  boundary  of  the 
model.  Almost  all  the  warming  at  the  earth's  surface  is  due  to  the 
absorption  of  incoming  solar  radiation.  Cooling  occurs  due  to  both 
the  net  loss  of  long-wave  radiation  and  the  loss  of  heat  through  evapor¬ 
ation,  the  latter  component  being  released  into  the  atmosphere  as  latent 
heat  of  condensation  when  precipitation  occurs.  The  earth's  surface 
may  also  be  cooled  or  heated  due  to  the  transfer  of  sensible  heat,  de¬ 
pending  on  whether  the  surface  is  warmer  or  cooler  than  the  adjacent 
atmosphere . 

The  zonally  averaged  January  heat  budgets  at  the  surface,  and  their 
components,  are  shown  in  Fig.  18  for  the  Mintz-Aru'cawa  model,  in  Fig.  19 
for  the  GFDL  model  discussed  by  Holloway  and  Manabe  (1971),  and  in  Fig. 
20  for  the  data  of  Budvko  (1963) .  Figures  18  and  19  include  not  only 
the  long-  and  short-wave  components  but  also  the  net  radiation  (the  sum 
of  the  two  parts),  which  is  the  only  radiative  quantity  available  from 
Budyko.  On  comparing  the  individual  components  of  the  heat  budget  in 
the  three  figures,  it  is  evident  that  the  Mintz-Arakawa  model  has  the 
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largest  components  of  both  downward  and  upward  heat  transfer.  The  net 
surface  flux  in  the  Mintz-Arakawa  model  is  larger  than  that  found  in 
either  of  the  ocher  studies  and  accordingly  yields  greater  radiative 
warming  at  the  earth's  surface.  In  addition,  whereas  in  the  other  two 
figures  there  is  a  loss  of  heat  at  the  surface  at  all  latitudes  due  to 
sensible  heat,  in  the  Mintz-Arakawa  model  at  certain  latitudes  some 
amount  of  warming  occurs.  These  apparent  overestimations  of  heating 
are  to  some  extent  balanced  by  the  fact  that  the  Mintz-Arakawa  model 
also  has  the  largest  values  of  surface  heat  loss  due  to  evaporation. 

The  net  surface  heat  budget  over  the  entire  globe  is  small  in  all  three 
cases,  but  the  zonal  structure  is  such  that  the  Mintz-Arakawa  model  has 
larger  negative  values  in  the  mid-northern  latitudes  and  larger  positive 
values  in  mid-southern  latitudes,  with  a  steeper  north-south  gradient. 

The  heat  balance  of  an  earth-atmosphere  column  as  represented  in 
Fig.  21a  can  be  expressed  (in  the  notation  of  Newell  et  al.,  1969)  as 

^EA  =  (E  "  P)  +  <S0  +  SL  +  SA)  +  Div  (h  +  V  (D 

where  RN£A  =  net  radiation  across  upper  surface 

E  =  heat  lost  by  ground  surface  due  to  latent  heat  of 
evaporation 

P  =  heat  gained  by  atmosphere  due  to  latent  heat  of 
condensation  (precipitation) 

S0.  =  heat  stored  in  ocean,  land,  and  atmosphere, 

respectively  (as  temperature  changes) 

Fa,  Fq  =  vertically  integrated  (two-dimensional)  horizontal 
heat  flux  in  the  atmosphere  and  ocean, 
respectively ,  in  cal  cm  ^  day  ^  (see  Fig.  21b) 

The  heat  balance  for  the  atmospheric  column  alone  is 

(RN£A  -  R  )  +  SH  +  P  =  SA  +  Div  Fa  (2) 

where  R  =  net  radiation  downward  across  the  lower  surface 
s 

SH  =  sensible  heat  gained  by  the  atmosphere  from  the 
lower  surface 
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(b)  Horizontal  heat  transport 


(a)  Column  heat  balance 


Fig.  21-Schematic  representation  of  the  heat  balance 
for  an  earth-atmosphere  column. 

In  the  Mintz-Arakawa  model  all  the  quantities  appearing  in  these  equa¬ 
tions  are  known  except"*^,  FQ,  and  SQ.  The  JUm  SQ  +  Div"^  can  be 
found  from  SQ  +  SL  +  Div  fJ*-  R8  -  SH  -  E,  with  SL  defined  to  be  equal 
to  zero.  Either  of  these  equations  can  then  be  solved  to  find  Div~F~. 

The  zonally  averaged  January  heat  budget  terms  for  the  atmospheric 
column  are  shown  in  Figs.  22,  23,  and  24  for  the  Mintz-Arakawa  model, 
the  GFDL  model,  and  a  composite  of  Moller's  (1951)  condensation  data, 
Budyko’s  (1963)  sensible  heat  data,  and  London's  (1957)  and  Sasamori 
et  al.'s  (1972)  radiation  data.  Rather  than  plot  the  radiation  terms 
as  given  in  Eq.  (2),  we  have  plotted  the  distributions  of  solar  radia¬ 
tion  absorbed  by  the  atmospheric  column  and  net  long-wave  radiation 
lost  by  the  atmospheric  column.  Also,  the  atmospheric  storage  and 
horizontal  energy  flux  terms ,  SA  and  Div  F^  have  been  combined  to  form 
a  single  term  called  the  total  or  imbalance  term.  This  was  necessary 
because  the  separate  terms  were  only  available  for  the  Mintz-Arakawa  * 


lally  averaged  heat  budget  of  an  atmospheric  column  for  the  Mi ntz-Arakawa  model 
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model.  Also,  this  makes  it  possible  to  compare  our  results  with  those 
of  Fig.  6  of  Newell  et  al.  (1969),  who  similarly  treats  the  two  terms 
together,  and  Fig.  35  of  Sellers  (1965) ,  who  considers  mean  annual 
values  and  assumes  that  the  storage  term  is  therefore  zero.  Values 
for  the  atmospheric  heat  storage,  S  ,  as  calculated  by  the  Mintz- 
Arakawa  model  are  given  in  the  Appendix,  where  it  can  be  seen  that  S 
comprises  only  a  small  fraction  of  the  imbalance  term,  the  remainder, 

Div  Fa,  representing  the  heat  which  must  be  transported  meridionally 
by  the  atmosphere  to  effect  a  balance. 

In  comparing  these  five  figures  (Figs.  22,  23,  24  plus  Fig.  6  of 
Newell  et  al.  (1969)  and  Fig.  35  of  Sellers  (1965))  several  features 
may  be  noted.  Perhaps  most  striking  is  the  similarity  of  the  shapes 
of  the  condensation  and  imbalance  curves  in  each  figure  (AC  in  Sellers' 
figure).  This  is,  of  course,  because  all  the  other  components  vary 
much  more  slowly  with  latitude,  but  it  helps  to  emphasize  the  importance 
of  latent  heat  in  the  energy-transfer  process. 

The  imbalance  or  total  term  in  Fig.  24  and  from  Newell  et  al.  (1969), 
both  representing  values  derived  from  a  mixture  of  observed  and  theoret¬ 
ical  data,  do  not  appear  to  be  realistic.  They  both  indicate  that  the 
imbalance  is  predominantly  negative,  which  is  physically  impossible. 

When  globally  averaged  the  imbalance  should  go  to  zero  except  for  the 
small  storage  term.  Newell  et  al.  have  suggested  that  possible  errors 
in  the  precipitation,  RN£A>  or  sensible  heat  could  be  involved.  The 
general  circulation  models,  presented  in  Figs.  22  and  23,  are,  on  the 
other  hand,  constructed  to  conserve  energy  and  therefore  should  not 
have  this  problem.  It  may  be  noted,  however,  that  the  Mintz-Arakawa 
model  does  not  have  an  exact  energy  balance  due  to  a  combination  of 
factors,  to  be  discussed  later. 

On  comparing  the  atmospheric  heat  budget  of  the  Mintz-Arakawa 
model  with  the  other  heat  budgets,  we  note,  as  before,  that  the  Mintz- 
Arakawa  values  for  most  of  the  components  are  the  largest,  both  positive 
and  negative,  implying  a  more  vigorous  general  circulation. 


Newell  et  al.  (1969)  speak  of  boundary  flux,  which  we  assume  is 
the  sensible  heat  flux. 
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Figures  25  and  26  show  the  various  components  making  up  the  heat 
budget  for  the  whole  earth-atmosphere  system  for  the  Mintz-Arakawa 
model  and  from  tables  of  data  from  Newell  et  al.  (1969)  for  December, 
January,  and  February.  In  the  Mintz-Arakawa  model,  with  its  fixed- 
temperature  infinite-heat-capacity  ocean,  it  is  not  possible  to  com¬ 
pute  Sq  or  Div  Fq  separately,  but  only  their  sum.  To  complete  Fig.  25, 
therefore,  we  have  adopted  the  values  of  Sq  given  by  Newell  et  al.  in 
order  to  calculate  Div  Fq  for  the  Mintz-Arakawa  model.  The  imbalance 
term  in  these  figures  is  Div  (F^  +  Fq) ,  the  divergence  of  heat  trans¬ 
ported  by  the  ocean  and  atmosphere.  Since,  as  already  noted,  the 
atmospheric  heat  storage  term,  S^,  is  quite  small,  it  has  not  been 
plotted  in  Figs.  25  and  26.  There  is  fairly  good  agreement  between  the 
figures.  The  larger  positive  values  of  the  imbalance  term  in  the  high 
southern  latitudes  in  the  Mintz-Arakawa  model  are  due  to  the  overly 
high  net  radiation  in  that  region,  as  already  noted.  The  larger  nega¬ 
tive  values  in  the  middle  northern  latitudes  are  due  to  the  differences 
in  the  precipitation  patterns. 

The  horizontal  flux  of  heat  in  the  ocean  and  atmosphere  was  cal¬ 
culated  by  numerically  integrating  the  zonal  averages  of  Div  (Fq  +  F^) 
with  the  appropriate  area  weighting.  This  was  done  for  the  Mintz- 
Arakawa  model  by  starting  at  the  south  pole  and  integrating  northward, 
in  the  manner  described  by  Newell  et  al.  (1969).  They,  however,  re¬ 
distributed  the  slight  residual  at  the  north  pole  among  the  largest 
flux  values,  wh51e  we  did  not.  The  result  for  the  Mintz-Arakawa  model 
and  Newell  et  al.  are  shown  in  Fig.  27.  Fairly  good  agreement  exists, 
with  the  Mintz-Arakawa  model  requiring  a  larger  horizontal  flux  in  the 
southern  hemisphere  and  a  considerably  smaller  flux  in  the  northern 
hemisphere  than  that  indicated  by  Newell  et  al. 

Finally,  we  will  consider  the  global  averages  of  the  energy  balance. 
From  Eqs.  (1)  and  (2),  we  have 


RNEA  =  (£  -  p)  +  Sq  +  5L  +  SA  +  Div  (PA  +  Fq) 

4  ■  260  -  258  -  1  +  0  +  5  +  (23  -  25) 


(1 1 ) 


ewall  et  al .  (1969). 
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with  the  numerical  values  of  the  global  averages  of  the  heat  budgets 
in  ly  day  .  Similarly, 


RNea  -  Rs  +  SH  +  p  =  SA  +  DivTA 
4  -  242  +  8  +  258  =  5  +  23 


(2') 


Equations  (1')  and  (2')  indicate  that  there  is  a  net  4  ly  day-1  enter¬ 
ing  the  atmosphere  across  the  upper  boundary  (RNgA) ,  and  24  ly  day-1 
across  the  lower  boundary  (-R,  +  SH  +  P).  Of  the  total  of  28  ly  day"1, 
only  5  ly  day  are  stored  in  the  atmosphere  as  heat.  The  other  23 
ly  day  are  unaccounted  for,  and  since  Eq.  (2')  forces  the  atmospheric 
system  to  be  in  balance,  these  23  ly  day-1  appear  in  the  horizontal 
atmospheric  heat-transfer  term,  which  should  vanish  globally.  At  the 
present  time,  however,  frictional  heating  of  the  earth's  surface  is  not 
included  in  the  model,  and  it  would  appear  that  if  this  surface  heating 
is  introduced,  part  of  the  missing  heat  will  be  accounted  for.  In  addi¬ 
tion,  the  use  of  single  precision  arithmetic  may  contribute  to  the 
imbalance.  Finally,  some  of  the  heat  energy  may  be  converted  into  other 
forms  of  energy.  The  net  global  horizontal  transport  of  heat  in  the 
oceans  also  should  be  zero.  Since  the  sum  of  the  terms  in  Eq.  (1'), 
excluding  the  divergence,  is  quite  small,  the  non-zero  value  found  for 
Div  Fa  forces  Div  FQ  to  be  nonzero  if  the  equation  is  to  be  balanced. 

It  may  be  noted  that  this  global  average  error  is  3  percent  of  the 
incident  solar  radiation. 

The  global  heat  budget  of  the  Mintz-Arakawa  model  for  January  is 
shown  in  Fig.  28  and  compared,  where  possible,  with  London's  (1957) 
values  for  the  mean  annual  case.  The  largest  differences  are  in  the 
nonradiative  heat  transfer  occurring  across  the  lower  surface.  The 
Mintz-Arakawa  model  has  35  units  of  latent  heat  (where  the  solar  ra¬ 
diation  incident  on  the  top  of  the  atmosphere  represents  100  units) , 
and  1  unit  of  sensible  heat,  while  London  has  18.5  units  of  latent 
heat  and  11  units  of  sensible  heat.  In  addition,  London  has  forced 
the  left-hand  side  of  Eq.  (2)  to  be  zero,  while  the  Mintz-Arakawa 
model  has  a  3  percent  imbalance  as  noted  above,  and  also  a  small 


a  k  _  . 


28-The  global  heat  budgets  from  the  Mintz-Arakawa  model  for  January  (large  numbers) 
and  from  London's  mean  annual  data  (small  numbers). 
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amount  of  atmospheric  heat  storage  in  the  form  of  increases  in  the 

evaporated  moisture  and  temperature.  Globally,  the  solar  radiation 

absorbed  by  the  earth's  surface  in  the  Mints-Arakawa  model  is  about 

5  units  larger  (52  units  versus  47  units)  than  in  London's  results 

The  long-wave  radiation  lost  to  space  is  approximately  the  same,  while 

that  lost  from  the  earth's  surface  is  slightly  larger  in  the  Mints- 
Arakawa  model. 


-52- 


VI.  CONCLUDING  REMARKS 


In  this  report  we  have  examined  the  radiation  and  heat-budget  por¬ 
tions  of  the  Mintz-Arakawa  model  in  order  to  identify  areas  not  in 
agreement  with  the  theoretical  and  observed  values  available  for  com¬ 
parison.  In  general  those  features  of  the  model  considered  here  appear 
to  be  fairly  realistic.  In  addition,  as  noted  in  Sec.  V,  some  of  the 
heat  budgets,  such  as  that  of  Newell  et  al.  (1969),  based  on  a  combina¬ 
tion  of  theoretical  and  observed  data,  produce  unrealistic  values  for 
the  horizontal  heat  flux.  Therefore,  perhaps  we  should  not  be  too 
concerned  that  the  various  budget  terms  derived  from  the  Mintz-Arakawa 
model  deviate  as  they  do  from  the  terms  in  the  other  budgets.  It  would 
also  appear  that  some  of  the  departures  can  be  explained  in  terms  of 
deficiencies  in  the  moisture  and  cloud  parameterizations  of  the  model, 
but  improvements  in  this  area  must  await  modifications  of  the  model's 
hydrologic  cycle. 

There  are,  however,  some  areas  where  the  model  can  be  immediately 
improved.  It  appears  that  a  slightly  lower  solar  constant  is  more 
realistic,  and  that  a  representation  of  ozone  absorption  is  desirable. 

In  addition,  dividing  the  solar  spectrum  at  a  wavelength  other  than 
0.9  y,  as  discussed  in  Sec.  Ill,  may  improve  the  atmospheric  absorption. 
Adoption  of  these  modifications  should  improve  the  model  as  far  as  the 
solar  radiation  is  concerned.  For  the  long-wave  radiation  of  the  model, 
the  adoption  of  the  more  accurate  fit  to  the  Yamamoto  radiation  chart 
as  developed  at  UCLA  should  prove  advantageous.  We  also  hope  to  ex 
amine  the  radiation  in  the  model  by  using  observed  climatological  param¬ 
eters  rather  than  those  generated  by  the  model,  in  order  to  determine 
whether  or  not  the  observed  discrepancies  are  due  to  the  hydrological 
parameters  used  in  the  model,  as  has  been  suggested  here. 

* 

M.  E.  Schlesinger,  private  communication. 
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Appendix 

NUMERICAL  VALUES  FROM  THE  mintz-arakawa  model 


Zonal  and  global  averages  of  the  quantities  below,  as  generated  by 

the  Mintz-Arakawa  model  in  a  January  simulation,  are  tabulated  on  the 
following  pages. 


Table  Notation 

Quantity 

Units 

CL 

Cloudlneaa 

fractional 

W.V. 

Precipitable  water  vapor 

-2 

gm  cm 

S.  Alb. 

Surface  albedo 

fractional 

T 

g 

Ground  temperature 

°C 

E 

Evaporation  from  aurface 

mm  day  ^ 

P 

Precipitation 

mm  day  ^ 

s 

0 

Solar  radiation,  top  of  atmosphere 

ly  day  * 

S  -  Aba.  S. 

0 

Solar  radiation  reaching  troposphere 

ly  day  1 

Aba.  S. 

Solar  radiation  abaorbed  in  atratoBphere 

ly  day-1 

Refl. 

Total  reflected  aolar  radiation 

ly  day  ^ 

Aba.  T. 

Solar  radiation  abaorbed  in  tropoBphere 

ly  day  * 

S4 

Solar  radiation  reaching  ground 

ly  day  ^ 

S 

g 

Solar  radiation  absorbed  by  ground 

ly  day 

P.  Alb. 

Planetary  albedo 

fractional 

Ro 

Net  long-wave  radiation,  tropopauae 
(poaitive  upward) 

ly  day  * 

R4 

Net  long-wave  radiation,  ground  surface 

ly  day  * 

(poaitive  upward) 

R0-R4 

Long-wave  flux  divergence 

ly  day  * 

^ea 

Net  radiation,  tropopauae 

ly  day  * 

(poaitive  downward) 

R 

8  | 

Net  radiation,  ground  (positive  downward) 

ly  day  ^ 

E(heat) 

Latent  heat  loss  from  evaporation 

ly  day  1 

P (heat) 

Latent  heat  gain  from  precipitation 

ly  day  1 

SH 

Sensible  heat  (poaitive  upward) 

ly  day 

SA  ^ 

Heat  storage,  atmosphere 

ly  day  1 

Div  Fa 

A  m 

Atmospheric  horizontal  heat-flux  divergence 

ly  day  1 

so  +  Dlv  Fo 

Heat  transferred  to  ocean 

ly  day  1 

s0 

Heat  storage,  ocean  (from  Newell  et  al.) 

ly  day 

Div  Fq 

Oceanic  horizontal  heat-flux  divergence 

ly  day  ^ 
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ZuNAL,  HEMISPHERIC  (NH,  SH) ,  AND  GLOBAL  AVERAGES  FROM  MINTZ-ARAKAWA  MODEL 


Lac. 

Rail. 

Aba.  T. 

S4 

S 

8 

90N 

0 

0 

0 

0 

86N 

0 

0 

0 

0 

82N 

0 

0 

0 

0 

78N 

0 

0 

0 

0 

74N 

0 

0 

0 

0 

70N 

0 

0 

0 

0 

66N 

10 

2 

8 

4 

62N 

32 

7 

28 

15 

58N 

59 

14 

48 

30 

54N 

88 

21 

70 

48 

SON 

121 

29 

92 

64 

46N 

149 

37 

120 

88 

42N 

162 

44 

150 

127 

38N 

177 

50 

189 

166 

34N 

183 

55 

242 

214 

30N 

181 

60 

306 

269 

26N 

192 

67 

352 

308 

22N 

206 

76 

388 

340 

1SN 

196 

85 

441 

395 

14N 

201 

93 

472 

431 

ION 

205 

98 

504 

470 

6N 

220 

104 

528 

493 

2N 

227 

111 

556 

520 

2S 

238 

117 

578 

540 

6S 

247 

121 

602 

562 

10S 

252 

125 

625 

583 

14S 

257 

127 

647 

601 

18S 

265 

130 

662 

613 

22S 

285 

134 

665 

606 

26S 

314 

136 

643 

589 

30S 

370 

138 

593 

540 

34S 

384 

137 

574 

533 

38S 

384 

134 

576 

537 

42S 

384 

132 

581 

537 

46S 

398 

133 

560 

516 

SOS 

416 

132 

530 

490 

54S 

413 

128 

531 

485 

58S 

433 

130 

494 

452 

62S 

416 

125 

529 

462 

66S 

403 

118 

569 

476 

70S 

548 

108 

644 

353 

74S 

566 

99 

672 

367 

78S 

752 

103 

734 

194 

82S 

780 

104 

733 

179 

86S 

822 

109 

701 

139 

90S 

822 

n: 

697 

138 

NH 

157 

59 

285 

215 

SH 

353 

127 

601 

528 

Global 

255 

93 

443 

392 

P.  Alb. 

Ro 

R4 

R„  -  8, 

0  4 

™EA 

0 

303 

108 

195 

-304 

0 

300 

110 

190 

-301 

0 

300 

112 

188 

-301 

0 

302 

113 

189 

-303 

0 

305 

112 

193 

-306 

.087 

307 

108 

198 

-307 

.617 

307 

101 

207 

-303 

.582 

317 

101 

217 

-296 

.563 

329 

101 

228 

-286 

.555 

339 

100 

239 

-271 

.562 

348 

104 

245 

-257 

.540 

365 

115 

249 

-242 

.483 

389 

132 

258 

-220 

.447 

416 

147 

269 

-202 

.401 

448 

159 

289 

-181 

.353 

486 

179 

307 

-158 

.337 

522 

190 

332 

-148 

.328 

544 

188 

356 

-127 

.288 

565 

192 

372 

-85 

.275 

574 

191 

383 

-50 

.264 

578 

190 

388 

-9 

.267 

579 

184 

395 

19 

.263 

580 

176 

404 

52 

.264 

577 

171 

406 

82 

.263 

575 

165 

409 

108 

.261 

575 

163 

413 

133 

.259 

573 

161 

411 

157 

.261 

570 

161 

408 

175 

.276 

568 

169 

399 

173 

.300 

556 

166 

389 

170 

.351 

532 

155 

377 

148 

.361 

509 

148 

362 

164 

.362 

486 

135 

351 

187 

.363 

465 

129 

336 

204 

.378 

439 

117 

322 

212 

.398 

414 

108 

306 

210 

.400 

396 

107 

289 

221 

.423 

375 

95 

281 

210 

.412 

372 

99 

273 

218 

.401 

371 

110 

261 

226 

.538 

372 

133 

239 

91 

.545 

370 

150 

221 

98 

.711 

341 

125 

216 

-42 

.728 

332 

119 

213 

-48 

.762 

319 

107 

213 

-70 

.760 

319 

105 

214 

-67 

.363 

464 

155 

309 

-150 

.345 

499 

144 

355 

158 

.354 

481 

150 

332 

4 
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ZONAL,  HEMISPHERIC  (NH,  SH) ,  AND  GLOBAL  AVERAGES  PROM  M1NTZ-ARAKAWA  MODEL 


uc. 

R 

a 

E  (hast) 

P(haat) 

SH 

SA 

d1vTa 

SQ  +  DivT0 

-so 

Div  *o 

90N 

-108 

-2 

7 

-75 

-2 

-261 

-31 

0 

-31 

86N 

-110 

6 

62 

-8A 

1 

-213 

-32 

9 

-23 

82N 

-112 

A 

A3 

-93 

2 

-238 

-25 

16 

-9 

78N 

-113 

13 

A3 

-83 

12 

-2A1 

-AA 

30 

-1A 

7  AN 

-112 

38 

56 

-AO 

11 

-188 

-111 

AO 

-71 

70N 

-108 

30 

61 

-6A 

2 

-202 

-7A 

A5 

-29 

66N 

-97 

20 

67 

-72 

-1 

-208 

-A6 

53 

7 

62N 

-86 

2A 

68 

-62 

6 

-210 

-A8 

59 

11 

5  BN 

-71 

A2 

80 

-30 

13 

-177 

-82 

67 

-15 

5AN 

-52 

3A 

103 

-A1 

1A 

-169 

-A5 

78 

33 

50N 

-AO 

68 

1A9 

-8 

17 

-91 

-99 

00 

-9 

A6N 

-27 

126 

215 

31 

19 

15 

-185 

105 

-80 

A2N 

-5 

22A 

261 

76 

20 

10A 

-305 

125 

-180 

3  BN 

19 

320 

30A 

106 

1A 

178 

-A08 

1A2 

-266 

3AN 

55 

331 

238 

111 

7 

108 

-387 

150 

-237 

30N 

90 

299 

132 

95 

1 

-20 

-30A 

1A5 

-159 

26N 

118 

332 

100 

68 

2 

-100 

-282 

128 

-15A 

22N 

152 

3A2 

126 

50 

2 

-106 

-2A0 

100 

-1A0 

1SN 

203 

396 

260 

27 

3 

-A 

-221 

80 

-1A1 

IAN 

2A0 

A15 

318 

19 

A 

A3 

-19A 

6A 

-130 

ION 

280 

A09 

A2A 

15 

A 

1AA 

-1A3 

AO 

-103 

6N 

309 

39  B 

A13 

13 

5 

130 

-102 

2A 

-78 

2N 

3AA 

370 

A59 

0 

A 

162 

-27 

8 

-19 

2S 

369 

350 

511 

-9 

A 

208 

29 

-A 

25 

6S 

397 

368 

531 

-21 

5 

216 

A9 

-15 

3A 

10S 

A20 

367 

A87 

-2A 

5 

170 

78 

-30 

A8 

IAS 

AAO 

380 

A72 

-3A 

5 

1A9 

93 

-52 

A1 

18S 

A52 

391 

AA7 

-3A 

5 

130 

95 

-80 

15 

22S 

A37 

359 

3A3 

-2 

A 

72 

80 

-108 

-28 

26S 

A23 

3A2 

300 

18 

3 

62 

63 

-122 

-59 

30S 

385 

319 

221 

62 

2 

A2 

A 

-135 

-131 

3AS 

385 

309 

17A 

63 

2 

9 

1A 

-1A0 

-126 

388 

A02 

287 

170 

25 

3 

-26 

90 

-1A0 

-50 

A2S 

AO  8 

210 

161 

1 

A 

-A6 

195 

-139 

56 

A6S 

399 

161 

178 

-18 

1 

-31 

256 

-132 

12A 

50S 

382 

11A 

187 

-23 

0 

-10 

290 

-125 

165 

5AS 

378 

75 

159 

-29 

-A 

-26 

332 

-105 

227 

58S 

357 

2A 

117 

-37 

-6 

-6A 

370 

-78 

292 

62S 

363 

10 

79 

-A5 

1 

-11A 

398 

-A8 

350 

66S 

366 

21 

52 

-A8 

3 

-1A2 

39A 

-27 

367 

70S 

220 

77 

106 

-39 

-1 

-63 

183 

-15 

168 

7AS 

217 

112 

129 

0 

-1 

8 

106 

-A 

102 

78S 

69 

A3 

52 

-22 

-2 

-81 

AS 

0 

A8 

82S 

60 

-'ll 

61 

-15 

-8 

-5A 

3A 

0 

3A 

86S 

32 

AA 

73 

-11 

-9 

-33 

0 

0 

0 

90S 

33 

3A 

60 

0 

-8 

-3A 

0 

0 

0 

NH 

100 

261 

227 

2A 

7 

-5 

-185 

82 

-103 

SH 

38A 

259 

288 

-7 

2 

51 

133 

-80 

53 

Global 

2A2 

260 

258 

8 

5 

23 

-26 

1 

-25 
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